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£j\ Abstract 

S: 

We demonstrate that perturbative QCD allows one to calculate the absolute 

cross section of diffractive exclusive production of photons at large Q 2 at 
"xl" ' 

;> ■ 

| HERA, while the aligned jet model allows one to estimate the cross section for 

in 

intermediate Q 2 ~ 2GeV 2 . Furthermore, we find that the imaginary part of 

O ' 



the amplitude for the production of real photons is larger than the imaginary 
part of the corresponding DIS amplitude by about a factor of 2, leading to 
the prediction of a significant counting rate for the current generation of 
experiments at HERA. We also find a large azimuthal angle asymmetry in ep 
scattering for HERA kinematics which allows one to directly measure the real 

X ' 

$—1 ' part of the DVCS amplitude and hence the nondiagonal parton distributions. 

PACS: 12.38.Bx, 13.85.Fb, 13.85.Ni. 
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I. INTRODUCTION 

Recent data from HERA has spurred great interest in exclusive or diffractive direct 
production of photons in e — p scattering (DVCS- deeply virtual Compton scattering) as 
another source to obtain more information about the gluon distribution inside the proton 
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for nonforward scattering. In recent years studies of diffractive vector meson production 
and deeply virtual Compton scattering has greatly increased our theoretical understanding 
about the gluon distribution in nonforward kinematics and how it compares to the gluon 
distribution in the forward direction. For a less than complete list of recent references see 
Ref. §0. 



Exclusive diffractive virtual Compton processes at large Q , first investigated in ||12 



offer a new and comparatively "clean" |T3| way of obtaining information about the gluons 



inside the proton in a nonforward kinematic situation. We are interested in the production 
of a real photon compared to the inclusive DIS cross section. The exclusive process is 
nonforward in its nature, since the photon initiating the process is virtual (q 2 < 0) and the 
final state photon is real, forcing a small but finite momentum transfer to the target proton 
i.e forcing a nonforward kinematic situation as we would like. We will show that pQCD can 
be applied to this type of exclusive process although we will not give a formal proof on the 



level comparable to the DIS case |14|]. Such a proof can be found in Ref. |p~5|| . 

The paper is organized in the following way. In Sec. [TJ we estimate the amplitude in the 
normalization point Ql ~ 2GeV 2 using the aligned jet model approximation and conclude 
that for such Q 2 the nondiagonal amplitude is larger than the diagonal one by a factor of 
~ 2. In Sec. |Tl| we calculate the imaginary part of the amplitude for 7* + p — > 7 + p in the 
leading order of the running coupling constant a s and compare it to the imaginary part of 
the amplitude in DIS in the same order. In Sec. [IV] we argue that at sufficiently small x the 
t-dependence of the cross section should reflect the interplay of hard and soft physics typical 
of diffractive phenomena in DIS. Namely, that at fixed x and increasing Q 2 , hard physics 
should tend to occupy the dominant part of the space of rapidities. In contrast to this, at 
fixed Q 2 and decreasing x, hard physics should occupy a finite range of rapidities which 
increases with Q 2 - ~ In -j^p with /3 ~ 0.1 — 0.2 at the HERA energy range due to the 
QCD evolution, and that soft QCD physics occupies the rest of the phase space. In Sec. [V] 
we give the total cross section of exclusive photon production and give numerical estimates 
of the DVCS production rate at HERA and find that such measurements are feasible for 



the current generation of experiments. We also show the feasibilty of directly measuring 
the real part of the DVCS amplitude and hence, at least, the shape of the nondiagonal 
parton distributions through a large azimuthal angle asymmetry in ep scattering for HERA 
kinematics. Sec. |V11| finally contains concluding remarks. 

II. THE AMPLITUDE FOR DIFFRACTIVE VIRTUAL COMPTON 
SCATTERING AT INTERMEDIATE Q 2 

Similar to the case of deep inelastic scattering, in real photon production it is possible 
to calculate within perturbative QCD the Q 2 evolution of the amplitude but not its value 
at the normalization point at Qq ~ few GeV 2 where it is given by nonperturbative effects. 



Hence we start by discussing expectations for this region. It was demonstrated in |17j that 
the aligned jet model |16| coupled with the idea of color screening provides a reasonable 
semiquantitative description of F 2 n{x < 10~ 2 , Qq). In this model the virtual photon interacts 
at intermediate Q 2 and small x via transitions to a qq pair with small transverse momenta 
- ^o,i (\^oi) ~ 0.15GeV 2 ) and average masses ~ Q 2 which thus carry asymmetric fractions 



of the virtual photon's longitudinal momentum. Due to large transverse color separation, 



b ~ 2v/2/3r 7r , the aligned jet model components of the photon wave function interact strongly 
with the target with the cross section a tot ( Ll AJM"—N) w a tot (7rN). Neglecting contributions 
of the components of the 7* wave function with smaller color separation, one can write 



&tot{l*N) using the Gribov dispersion representation [18[] as |T7 



a tot {l N) = - (Qi + M 2 ) 2 dM ' (1) 

where the factor M 2 in the nominator is due to the overall phase volume, R e+e (M 2 ) = 
^(e+e^^+zT- ^ • ^ ne f ac t° r 3 ls ^ e f rac ti° n °f the whole phase volume occupied by the 
aligned jet model , and the factor 1/{Q 2 + M 2 ) 2 is due to the propagators of the photon in 
the hadronic intermediate state with mass 2 equal M 2 . Based on the logic of local quark- 
hadron duality (see e.g. [113 and references therein) we take the lower limit of integration 
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Mq ~ 0.5GeV 2 < m 2 . In the case of real photon production the imaginary part of the 
amplitude for t = is 

1 W + S^ + ^ - £ Q °^'^^^ M, (2 ) 

with s = 2qorriN being the flux factor. The only difference from Eq. [I] for <J to t{l* + N) is the 
change of one of the propagators from 1/(Q 2 + M 2 ) to 1/M 2 - here go is the energy of the 
virtual photon in the rest frame of the target. 

Approximating R e+e (M 2 ) as a constant for the Q 2 range in question (we understand 
this in the sense of a local duality of the hadron spectrum and the qq loop) we find 



_ ImA(Y + N-^r + N) t =o Q 2 _ 1( n2/ 2 ( . 

R = Tzrn-L i at aa — = ^2 I ^2 ln (l + QMo). (3) 



imA(7* + -»• 7 + iV) 4=0 Q 2 + M, 



In the following analysis we will take Qq for the perturbative QCD evolution as 2.6 GeV 
to avoid ambiguities. It is easy convince oneself that for M5 2 ~ 0.4 ~ 0.6 GeV 2 and Q 2 ~ 
2 — 3 GeV 2 Eq.|] leads to R ~ 0.5. A similar value of i? has been obtained within the 
generalized vector dominance model |H| As we will see below QCD evolution leads to a 
strong increase of Q 2 ImA{ r y* + iV — > 7 + N) t=0 with increase of Q 2 for fixed x. However it 
does not change appreciably the value of R. 

III. THE AMPLITUDE FOR EXCLUSIVE REAL PHOTON PRODUCTION AT 

LARGE Q 2 . 

The process of exclusive direct production of photons in first nontrivial order of a s ln %i 
at small Xsj can be calculated (see Fig. P as the sum of a hard contribution calculated within 



the framework of QCD evolution equations |25] and a soft contribution which we evaluated 



above within the aligned jet model. The hard contribution can be described through a two 
gluon exchange of a box diagram with the target proton. In order to calculate the imaginary 
part of the amplitude, we need to calculate the hard amplitude from the box as well as the 
gluon-nucleon scattering plus the soft aligned jet model contribution. Let us first give a 
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FIG. 1. Leading contribution to DVCS at small x. 

general expression for the imaginary part of the amplitude and then proceed to deal with 
the gluon-nucleon scattering, followed by the calculation of the box diagrams. 

First, let us discuss the hard contribution which actually dominates in the considered 
process. To account for the gluon-nucleon scattering, we work with Sudakov variables for the 
gluons with momenta p\ and P2 attaching the box to the target and the following kinematics 
for the gluon-nucleon scattering: 

Pi = aq' + Xip' + p t , d 4 pi = -dadxid 2 p t , (4) 

where q' and p' are light-like momenta related to p, q the momenta of the target proton and 
the probing virtual photon respectively, by: 

q = q' - xp, p = p' + ^—q', 

2p'q' 

s = 2pq = 2p'q' — xp 2 , (5) 

with x being the Bjorken x and x\ the proton momentum fraction carried by the outgoing 
gluon. Equivalent equations to Eq. |4] apply for pi with the only difference being that x\ 
is replaced by X2, the momentum fraction of the incoming gluon, signaling that there is 
only a difference in the longitudinal momenta but not in the transverse momenta. This 
fact will shortly become important. Furthermore there is a simple relationship between x\ 
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and X2 '■ A = xi — X2 = const, following from the kinematics of the considered reaction 



2l| where A is the asymmetry parameter or skewedness of the process under consideration. 
Therefore one is left with just the integration over p\ since p 2 cannot vary independently 
of p\. Being exclusively interested in the small x region, one can safely make the following 
approximations: s = 2pq ~ 2p'q' and p' ~ p. Since we are working in the leading a s In Q 2 
approximation, neglecting corrections of order a s , the main contribution comes from the 
region p\ << Q 2 , hence the contribution to the imaginary part of the amplitude simplifies 
considerably. First, since \p 2 \ = \ax±s + p 2 \ << Q 2 , one has a << 1 and the polarization 
tensor of the propagator of the exchanged gluon in the light-cone gauge q'^A^ = becomes, 



see Ref. [18 



4a * (6) 
p'q' v ; 

In other words it is enough to take the longitudinal polarizations of the exchanged gluons 
into account. 

Using Eq. |^ one obtains the following expression for the total contribution of the box 
diagram and its permutations: 

ImA = f ^^ImAfPlmAfPd.Md^), (7) 

where lmA a ^ p ^ = ImA°^,(j* g — > qq) is the sum of the box diagrams, ImA ^^ is the 
amplitude for the gluon-nucleon scattering, a,b are the color indices and the overall tensor 
structure has been neglected for now. The usage of the imaginary part of the scattering 
amplitude and in particular limiting ourselves to the s-channel contribution as the dominant 
part in both the forward and the nonforward case (Eq. [7]) is correct (see Ref. for more 
details) as long as we restrict ourselves to the DGLAP region of small x and thus small t, 
where t = (pi — p 2 ) 2 is the square of the momentum transfered to the target. The real part 
of the amplitude will be evaluated below by applying a dispersion relation over the center 
of mass energy s. Using Eq. |5] and the Ward identity which is the same as in the Abelian 
case since the box contains no gluons i.e is color neutral: 
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<J P W = 0,< (P W = 0, (8) 

yielding 

ImAf^p'^ = ImA a ^p t ^p tv 
4(pg) 2 X1X2S 2 

one can rewrite Eq. |7| as: 

ImA f 1 dx x , 2 2 2 f sdad 2 p t 2 4ImA a x { pq x q a 

= / E{x x u A x u Q ,p t ,Q ) / 2 2 #£ Q *f , 10 

where we have used < pt^Ptu >= — ^dtivPt (the average over the transverse gluon polariza- 
tion) and defined the imaginary part of the hard scattering to be given by: 

1 ImA ab ^ 

Eix/x^A/x^Q^plQl) = --g* ^—5 abl (11) 

2 p X2S 

where the sum over repeated indices is implied. Up to this point we have just rewritten the 
equation for the imaginary part of the total amplitude but have not identified the different 
perturbative and non-perturbative pieces. In the case of a virtual photon with longitudinal 
polarization, this would be an easy task since the qq pair would only have a small space- 
time separation and we could follow the argument in Ref. |l|J^,|l^] stating that the box is 
entirely dominated by the hard scale Q and thus can unambiguously be calculated in pQCD. 
However, in our case we are dealing with a virtual photon which is transversely polarized 
and thus one can have large transverse space separations between q and q. The resolution to 
this problem can be found in the following way: one accepts that one has a contribution from 
a soft, aligned-jet-model-type, configuration and that there is no unambiguous separation 
of the amplitude in a perturbative and non-perturbative part up to a certain scale Qq. 
However, in the integration over transverse gluon momenta, one will reach a scale at which 
a clear separation into perturbative and non-perturbative part can be made and hence we 
can unambiguously calculate albeit not the imaginary part of the amplitude of the upper 
box but its InQ 2 derivative i.e its kernel convoluted with a parton distribution. At this 
point then, one can include the non-perturbative contribution of the aligned jet model into 
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the initial distribution of the imaginary part of the total amplitude and solve the differential 
equation in Q 2 . One obtains the following solution for the imaginary part fl25fl: 



ImA{x, Q 2 , Ql) = ImA(x, Q 2 ) + / \ -^P qg ( x /x 1 , A/x 1 )g(x 1 , x 2 , Q' 2 ), (12) 

JQq V Jx X\ 

where P qg is the evolution kernel |26| and starting from Ql the gluon distribution can be 



defined from Feynman diagrams in the leading a s In Q 2 approximation by realizing that in 
Eq. |l^ one can replace pi and p\ by p 2 t and one finds: 

f sdad 2 p t 4ImA a x ( pq x q a 2 

J {2tiYp 2 a 7 2 = g(xi,X2,Q), (13) 

where g is the nondiagonal parton distribution in general. Comparison of Eq. |10| with the 
Q CD-improved parton model expression for the total cross section of charm production given 
p4| shows that g in the case A = is the conventional, diagonal gluon distribution. 



in 



Note that the parton distribution which serves as an input in Eq. |T| has to be evolved 
over the Q 2 " ran g e covered by the Q' 2 integral which complicates the calculation. We will 
explain below how to deal with this issue in practical situations. 

At this point we would like to comment on equivalent definitions of nondiagonal parton 
distributions in the literature which differ by kinematic factors (see for example [|1|||7|). Eq. 
[131 corresponds to the definition used in 0] , however since it is given on the level of Feynman 
diagrams there are no ambiguities such as renormalization of bilocal operators and hence it 
provides an unambiguous definition of a nondiagonal parton distribution. 

For the non-perturbative input, ImA(x, Ql) we will be able to use the aligned jet 
model analysis of Sec. [Ij and the standard relation between ImA' y * p ^' 1 * p (x,Q 2 ,t = 0) and 
F 2p (x,Q 2 ): 

ImA^ p ^ tp (x, Q 2 , t = 0) = F2p ^ Q2) . (14) 

Following the discussion above, we now only need to calculate P qg to leading logarith- 
mic accuracy, in order to make predictions for the imaginary part of the whole amplitude. 
Therefore, let us now consider the box diagram where the two horizontal quark propagators 
are cut, corresponding to the DGLAP region i.e neglecting the u-channel contribution. 
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FIG. 2. Cut box diagram giving the kernel for the imaginary part of the DVCS amplitude. 

The kinematics (see Fig. ^ for the calculation of the cut box diagram, using Sudakov 
variables, is the following. The quark-loop momentum k is given by: 



k = acq' + j3p' + p t , d 4 k = -dadf3d 2 k tl 
where q' and p' are light-like momenta related to p, q by: 



p 

q = q' - xp, p x = p + ^y-jg', 
s = 2piq = 2p'q' — xp\. 



(15) 



(16) 



The momenta of the exchanged gluons, in light cone coordinates, are given by: 



Pl = ( Xl p + ,0,p t ), p 2 = (x 2 p + ,0,pt), 



(17) 



where we have assumed the transverse momentum of the proton to be zero. The probing 
transverse photon and the produced photon have the following momenta, again in light cone 
coordinates: 



<l = ( -•'■/'+, 77^—- t ), qi = (~~ 0, -5— , ()/). 



(18) 



2xp + 2xp^ 

P qg is calculated in the light cone gauge yielding the following result for the most general 
case lETI: 
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P qg (x/ Xl , A/xO = 4n*a^N F x{x f } + ( *\ 2 ^ (19) 

7T — A)^ 

The DIS kernel is analogous to Eq. [TJ3| except that A = and the kernel for real photon 
production is obtained for A = x. 

We now can proceed to calculate the total imaginary part of the amplitude from Eq. [L2| 
where we parameterize the gluon distribution at small x as: 

g(x u x 2 ,Q 2 )=A (Q 2 )xf l{Q2) . (20) 



We neglect the x 2 dependence for the moment ||28|| . The above parameterization is taken 
from CTEQ3L as well as the parameterization of a in terms of Q 2 in leading order: 

MQ 2 ) = exphO.7631 - 0.7241 In ( - 1.171n 2 ( 

\\n(Q /A) Un(Q /A)/ 



+0.534 In 3 



( X 

\\n(Q /A)) 1 



M Q 2 ) = -0.3573 + 0.3469 In ( £}9!fjL) - 0.3396 In 2 ( pj&^L) 

+0.09188 In 3 f i^^ V (21) 



Vln(Q /A)y 



with A, Q and a s given by: 



A = 0.177 GeV Qo = 1.6GeV a s = ^^Jj^j' (») 

where we have taken Nq = 3 and Np = 3. 

The ratio R of the imaginary parts of the amplitudes [|2?J is given by: 

R = ImA^f +p-^ 7* + p) 
JmA(7* + p — > 7 + p) 

We give in the x range from 10~ 4 to 10~ 2 and for a Q 2 of 3.5, 12 and 45 GeV 2 since this 
kinematic range is relevant at HERA. One might ask what about the contributions due to 



quarks. The answer is that the corrections are small |30| but for completeness we include 
them here. Eq. [12] is then augumented with a similar expression for the quark contribution 
where the kernel is now that of quark-quark splitting and the nondiagonal parton distribution 
is that of the quark: 
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ImA(x,Q 2 ,Q 2 ) = ImA(x,Q 2 ) + / ^_ / — [P ?g (x/xi, A/x 1 )g(x i , x 2 , Q 12 ) 

JQq V •'s 2-1 

+P 99 (x/xi, A/xMxx, x 2 , Q'% (24) 

where the general expression for the kernel, after a similar calculation as before, is found to 
be: 

x/xi — x 3 /xl — A/xi (x/xi + a; 2 /a; 2 )) 



Pqq(x/Xx, A/X\) = A.TX 2 OL — C p 



(25) 



xx{l - Ajxx){l - x/xx)+ 
and the + - prescription is the one used in Ref. 0. The quark distribution itself is also 
taken from CTEQ3L [[3TJ and given by: 



q(x 1 ,x 2l Q 2 ) = A (Q 2 )xf\ (26) 



with 



A (Q 2 ) = exp[0.1907 + 0.04205 In ( Pj^r ) + 0.2752 In 2 ( HQ,K) ^ 



Vln(Q /A)y ' \ln(Q /A)J 



-0.3171 In 



An(Q/A) V 



Vln(Q /A); J 

Ax = 0.465. (27) 

We chose Ax to be constant since it varies only between 0.4611 and 0.468 in the Q 2 range of 
interest, i.e. the error we make is almost negligible since the quark distribution themselves 
are small in the x-range considered. Furthermore, according to our discussion in Sec. ||, we 
chose the initial distribution for the imaginary part of the DVCS amplitude to be twice that 
of the initial distribution for the imaginary part of the DIS amplitude. In the evolved QCD 
part, the nonforward kinamtics are taken into account in the kernels of the QCD evolution 
equation, also the different Q 2 evolution of nondiagonal as compared to diagonal distribution 
has been taken into account as explained below. 

As the calculation with MATHEMATICA showed, the amplitude of the production of 
real photons is larger than the DIS amplitude over the whole range of small x and R turns 
out to be between 0.551, 0.573 and 0.57 for x = 10~ 4 , 0.541, 0.562 and 0.557 for x = 10~ 3 
and 0.518, 0.519 and 0.505 for x = 10~ 2 in the given Q 2 range. It has to be pointed out 
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that for a given Q 2 , the ratio is basically constant. Of course, the ratio R will approach 1/2 
as Q 2 is decreased to the nonperturbative scale since this is our aligned jet model estimate 
The reason for the deviation from R— 1/2 is due to the difference in the evolution kernels. 

It is worth noting that in the kinematics we discuss, the ratio is still rather sensitive to the 
nonperturbative boundary condition. For example, assuming the same boundary conditions 
for DVCS and DIS, would result in a reduction of R of about 20(10)% at Q 2 ~ 12(40) GeV 2 
and x ~ 10~ 3 

In Eq. [12| the median point of the integral corresponds to x\/2 ~ £2 ~ x. This is due to 
the mass of the qq in the quark loop being oc Q 2 . For such & Xi/x 2 the ratio of nondiagonal 
and diagonal gluon densities weakly depends on x 2 - Hence with an accuracy of a few percent 
we can approximate this ratio by its value at x\jx% = 2. Therefore, in the calculation of 



R, we used Eq. [20] and ^6] for both the diagonal and nondiagonal case but then multiplied 
the real photon result of the amplitude by a function f(Q 2 ) for each x and Q 2 to take into 
account the different evolution of the nondiagonal distribution as compared to the diagonal 



rQ 2 HO' 2 r 1 r/r, 

ImA(x,Q 2 ,Q 2 ) = ImA(x,Q 2 )+ ~^f(Q' 2 ) —P qg (x/x 1 ,A/x 1 )g(x 1 ,x 2 ,Q' 2 ). (28) 

JQg v Jx X\ 

The function was determined by using our modified version of the CTEQ-package and, start- 
ing from the same initial distribution and evolving the diagonal and nondiagonal distribution 
to a certain Q 2 . We thenO compared the two distributions at the value x 2 = Xi/2 = x for 
different x and then interpolated for the different ratios of the distribution in Q 2 for given 
x. For this median point the difference between the diagonal and nondiagonal gluon distri- 
bution is between 8 — 25% depending on the x and Q 2 involved and — 5% for the quarks 
(see the figures in Ref. P,[TT]j for more details). 

As far as the complete amplitude at small x is concerned, we can reconstruct the real 
part via dispersion relations p2]j23| , which to a very good approximation gives: 

ReA it d \n(xImA) 



1mA 2 d In ± 

X 
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(29) 



Meaning that since ImA can be fitted as x _1_<5 , r\ m |5 is independent of x to a good 
precision. Therefore, our claims for the imaginary part of the amplitude also hold for the 
whole amplitude at small x. This is due to the fact that within the dispersion representation 
of the amplitude over x the contribution of the subtraction constant becomes negligible at 
sufficiently small x. 

One also has to note that there is a potential pitfall since the QED bremsstrahlung - the 
Bethe-Heitler process, where the electron interacts with a proton via a soft Coulomb photon 
exchange and the real photon is radiated off the electron, can be a considerable background. 
As was shown by Ji |J, the Bethe-Heitler process will give a strong background at small t 
and medium Q 2 and x > 0.1. We will discuss this subject in more detail later on. 



IV. THE i-SLOPE OF THE j*N -> jN CROSS SECTION 

The slope of the differential cross section of the virtual Compton scattering - d oc 

exp(.Bt) is determined by three effects: (i) the average transverse size of the qq component of 
the 7* and 7 wave functions involved in the transition, (ii) the pomeron-nucleon form factor 
in the nucleon vertex, and (iii) Gribov diffusion in the soft part of the ladder. This leads to 
several qualitative phenomena. In the normalization point, qq configurations of an average 
transverse size, comparable to that of the p-meson, give the dominant contribution to the 
scattering amplitude, leading to a slope similar to that of the processes 7+p — > p,u+p. The 
contribution of the higher mass qq components is known to result in an enhancement of the 
differential cross section of the Compton scattering at t = by a factor « 2 as compared to 
the prediction of the vector meson dominance model with p,u,<j), J/ip intermediate states, 
see e.g. [j3^] . Since the diffraction of a photon to masses Mx > 1.3GeV has a smaller t slope 



than for transitions to p and uj, one could expect that the high mass contribution would lead 
to a t-slope of the Compton cross section being somewhat smaller than for the production 
of p, a;-mesons. However direct experimental comparison [[!!J of the slopes of the Compton 
scattering and the cu-meson photoproduction at (E^ nc ) ~ lOOGeV^ finds these slopes to be 
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the same within the experimental errors. Using these data, we can estimate the slope of 
the amplitude for diffractive photon production in DIS at HERA energies but at moderate 
Q-i.e. in the normalization point as 

B(s,Ql) = B Gomp .Scau.(so) + 2a'ln(-), (30) 

so 



where a' = 0.25GeV~ 2 , s = 200GeV 2 , and B Comp _ Scatt _{s Q ) = 6.9 ± 0.3GeV~ 2 |3|. Hence 
for HERA energies B(W = 200GeV, Q 2 ) ~ lOGeV" 2 . 

In another limit of large Q 2 and large enough x, say x ~ 10~ 2 , the dominant qq configu- 
rations have small a transverse size and the upper vertex does not contribute to the slope. 
Furthermore, the perturbative contribution occupies most of the rapidity interval and leaves 
no phase space for the soft Gribov diffusion. In this case, the slope is given by the square 
of the two-gluon form factor of the nucleon which corresponds to B — B gg ^ «4-f- 5GeV~ 2 



An interesting situation emerges in the limit of large but fixed Q 2 when the energy starts 

to increase. In this case, the perturbative part of the ladder has the length ~ ln(^^). 

Here k = x/xq, where Xq is the fraction x of the parent parton at a soft scale. For HERA 

kinematics k ~ 0.1 — 0.3 for Q 2 ~ 10 — 20GeV 2 and increasing with increasing Q 2 . This 

is consistent with the observation of an approximate factorization for diffraction in the case 

of high masses (M 2 > lOOGeV 2 , M 2 ^> Q 2 ) in the scattering of real and virtual photons 

observed at HERA p4] , namely 

1 dajr/N -> XN)(W, M x ) _ 1 daj^N -> XN)(W, My) 
atotilN) dtdM 2 ~ cxtoti^N) dtdM 2 ' ( } 

The observed slope for these processes is B ~ 7GeV~ 2 which is consistent with the presence 

of a cone shrinkage at the rate ~ 2a' \n(W 2 / M 2 ) as compared to the data at lower energies 

where smaller values of W 2 /M 2 were probed. Similarly we can expect that for virtual 

Compton scattering at large Q 2 , the slope will increase with decrease of x, at very small x, 

approximately as 

W 2 K W 2 W 2 k W 2 

B(W 2 , Q% w = B ggN + 2a\M—) - M^)0(— - -|), (32) 
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where Wq = 200GeV 2 . We take into account here that B gg ^ was determined experimentally 
from the processes at W 2 ~ W 2 . 

V. THE RATE OF EXCLUSIVE PHOTON PRODUCTION AT HERA 

To check the feasibility of measuring a DVCS signal against the DIS background, we will 
be interested in the fractional number of DIS events to diffractive exclusive photoproduction 
events at HERA in DIS given by: 

# 7 = TV^ T+ *=° x ~5 a W P), 33 

with 

| (7 . + ^ 7 + p) .f|g) (1 + , ;V ', (34) 
from applying the optical theorem and where R is the ratio of the amplitudes given by Eq. 



2 2 2 2 

23], and 77 = ReA/Im A is given by Eq. |29], t = — Pt — ~Pt with t min = — m ^ x ~ 0. 



Note, that even though we used the total DVCS cross section in Eq. |33] we only need the 
ratio of the imaginary parts of the amplitudes to calculate R 1 . A complete expression for 
DVCS will be given in the next section. Note that only da/dt(t = 0) is calculable in QCD. 
The t dependence is taken from data fits to hard diffractive processes. 
Using the fact that F 2 (x, Q 2 ) ~ ^i°ih^}Q one can rewrite Eq. |33|: 

^"4^B F2(x ' g2)(1+r/2) - (35) 

where t] 2 ~ 0.09 — 0.27 for the given Q 2 range. We computed _R 7 , the fractional number of 
events given by Eq. p5L for x between 10 -4 and 10~ 2 and for a Q 2 of 2, 3.5, 12 and 45 GeV 2 



where the numbers for F2 were taken from p5[. Based on our analysis of the previous section 



we use Eq. [JO] for Q 2 = 2 GeV 2 , assuming that for Q 2 = 3.5 GeV~ 2 the slope drops by about 
1 -T- 2 units as compared to Eq. ^ to account for the decrease of the transverse size of the 
gg-pair; for larger Q 2 we use Eq. |3T]. 
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We find R-y ~ 1.1 x 10~ 3 , 9.9 x 10~ 4 at x = 10" 4 , 10" 3 and Q 2 = 2GeV 2 ; # 7 ~ 1.07 x 
10" 3 , 9.3 x 1(T 4 at x = 1(T 4 , 10" 3 and Q 2 = 3.5GeV 2 ; Ry ~ 4.5 x 1(T 4 , 3.78 x 10~ 4 2.5 x 1(T 4 
at x = 10~ 4 , 10~ 3 , 1(T 2 and Q 2 = 12GeV 2 ; and finally i? 7 ~ 1.49 x 10" 4 , 1.04 x 10" 4 at 
x = 10 -3 , 10~ 2 and Q 2 = 45GeV 2 . As is to be expected, the number of events rises at small 
x since the differential cross section is proportional to the square of the gluon distribution 
and the total cross section is just proportional to the gluon distribution i.e. the ratio in Eq. 
|33| is expected to be proportional to the gluon distribution and this assumption is born out 
by our calculation and falls with increasing Q 2 since F 2 does not grow as fast with energy. 



VI. THE COMPLETE CROSS SECTION OF EXCLUSIVE PHOTON 

PRODUCTION 

In order to study whether the Bethe-Heitler or DVCS Process will be dominant in real 
photon production we need the expressions for the differential cross sections first. 

We find that the differential cross section for DVCS can be simply expressed through 
the DIS differential cross section by multiplying the DIS differential cross section by R 1 (see 
Eq. |35|) which was calculated in the previous section. One can see this by observing how F 2 
is related to cr t ot(l*p) as given in Sec. [V] and (J to t(l*p) to ctdvcs via -R 7 in the same section. 
We then find using Eq. |35| for i? 7 



dcr DV cs na s s 



'l + il-yfy-^Fi^Q^l+r?) (36) 



dxdyd\t\d(j) r 4R 2 Q 6 

with &dvcs — daD dt CS l f =o x j? using the same exponential t dependence as in the previous 
section and R being the ratio of the imaginary parts of the DIS to DVCS amplitudes as 
computed earlier. 

In writing Eq. [36| we neglected Fl(x,Q 2 ) - the experimentally observed conservation of 
s channel helicities justifies this approximation - and assumed F 2 ~ 2xF%. y = 1 — E'/E 
where E' ist the energy of the electron in the final state and <p r = <pN + 4>e, where 0at is the 
azimuthal angle between the plane defined by 7* and the final state proton and the x — z 
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FIG. 3. The azimuthal final proton and electron angle in the transverse scattering plane. 



plane and e is the azimuthal angle between the plane defined by the initial and final state 
electron and the x — z plane (see Fig. 

In case of the Bethe-Heitler process, we find the differential cross section at small t to be 



da B H 



a 3 sy 2 (l + (1 - yY 



x 



G\{t) + tGIM 

1 + T 



(37) 



dxdyd\t\dcj) r 7rQ 4 |*|(l — y) 
with t = \t\/Am 2 N , s being the invariant energy and y the fraction of the scattered elec- 
tron/positron energy. GEif) and GmW are the electric and nucleon form factors and we 
describe them using the dipole fit 

1*1 , 



G E (t) ^ G D (t) 



0.71' 



and G M (t) = fi p G D (t), 



(38) 



where fi p = 2.7 is the proton magnetic moment. We make the standard assumption that 
the spin flip term is small in the strong amplitude for small t. 

In order to write down the complete total cross section of exclusive photon production 
we need the interference term between DVCS and Bethe-Heitler. Note that in the case of 
the interference term one does not have a spinflip in the Bethe-Heitler amplitude, i.e. , one 
only has Fi(t), as compared to Eq. |3"7j containing a spinflip part, i.e. , ^(t). The appropriate 
combination of G^(t) and Gm^) which yields Fi(t) is 



G E (t) + SrG M {t) 



1 + 



1*1 
Ami. 



(39) 
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We then find for the interference term of the differential cross section, where we already 
use Eq. [!7|, 

d<vcs+BH _ ^W 3 sy(l + (1 - y^cos^e'^F^x, Q 2 ) 



(40) 



dxdyd\t\d<p r 2Q 5 ^f(\t\)^(l - y)R 

~G E (t) + rG M (ty 



x 



1 + r 

with the + sign corresponding to electron scattering off a proton and the - sign corresponding 
to the positron. The total cross section is then just the sum of Eq. PB| , |57| and 

A. dependence of Bethe-Heitler as compared to DVCS for different Q 2 

At this point it is important to determine how large the Bethe-Heitler background is 
as compared to DVCS for HERA kinematics, hence, in the following discussion, we will 
estimate the ratio D allowing a background comparison: 

_ < d(T DVC s+BH/dxdydt > _ 
< ddBH / dxdydt > 



with < ... >= J 2?r d(j) T . Using the expressions from Sec. [VI], we compute D and find D > 1 
[see Figs. f|a and|6|a] for relatively small y and 0.1 < t < 0.6 with the given values of x and Q 2 
considered. Note, however, that this does not mean that the case for DVCS is hopeless. As 
it turns out, it is rather advantageous to have D < 1 when looking at the interference term 
which we will do next. Also note that there is a very strong energy dependence of D which 
extends the range where DVCS is significantly larger than Bethe-Heitler with increasing 
energy. 

It is convenient to illustrate the magnitude of the intereference term in the total cross 
section by considering the asymmetry for proton and either electron or positron to be in the 
same and opposite hemispheres ( we omit the rather cumbersome explicit expression but 



the reader can easily deduce it from Eq. p6,37 and 



^ _ J-tt/2 d(f) r da DV cs+BH - fnj 2 2 d<f) r da DV cs+BH 

Jo™ d(f) r d<JDVCS+BH 



in other words, one is counting the number of events in the upper hemisphere of the detector 
minus the number of events in the lower half, normalized to the total cross section. Fig. 
|5]a,b and [7]a,b show A for the same kinematics as above and we find that the asymmetry is 
fairly sizeable already for small t and is strongly dependent on the energy. Due to this fairly 
large asymmetry, one has a first chance to access nondiagonal parton distributions through 
this asymmetry. We will discuss A in more detail, in particular its energy dependence, in 
an upcoming paper. 

Note, there is an increased experimental difficulty to measure DVCS if the recoil proton 
is not detected in other words if t is not directly measured. However there is a simple, 
practical way around this problem which we will discuss next. 



B. DVCS alternative to tagged proton in the final state 

Another interesting process, which can be studied in the context of DVCS, is the one 
where the nucleon dissociates into mass "X" - 7* + p — > 7 + X. Perturbative QCD is 
applicable in this case as well. In particular the following factorization relation should be 
valid at sufficiently large Q 2 : 

The big advantage of the dissociation process as compared to the process where the target 
proton stays intact is that the Bethe-Heitler process is strongly suppressed for inelastic 
diffraction at small t due to the conservation of the electro-magnetic current, hence the 
amplitude is multiplied by an additional factor ^J~\t\ which is basically for the Bethe- 
Heitler process. Thus, the masking of the strong amplitude of photoproduction is small in 
this case. Since there is already data available on J/ip production, this quantity can give 
us information on how different the slopes for the production of massless to massive vector 
particles are, providing us with more understanding on how different or similar the exact 
production mechanisms are. Note that the ratio of the total dissociative to elastic cross 
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section of p meson production is found to be about 0.65 at large Q 2 which is basically 
of 0(1). The same should hold true for J/ip production and in fact this ratio should be 
a universal quantity. This is due to the fact that one has complete factorization, hence 
the hard part plus vector meson is essentially a point and thus for the soft part, is does 
not matter what kind of vector particle is produced. The above said implies for Eq. 
that it also should be of order unity, implying that the order of magnitude of the fractional 
number of events for real photon production to DIS remains unchanged even though the 
actual number of _R 7 might decrease by as much as 35%. 



VII. CONCLUSIONS 

In the above said we have shown that pQCD is applicable to exclusive photoproduction 
by showing that the ratio of the imaginary parts of the amplitudes of DIS to a real photon 
is calculable in pQCD after specifying initial conditions since the derivative in energy of 
the hard scattering amplitudes can be unambiguously calculated in pQCD and all the non- 
perturbative physics can then be absorbed into a parton distribution. We wrote down an 
evolution equation for the imaginary part of the amplitude, which can be generalized to the 
complete amplitude at small x, and solved for the imaginary part of the amplitude. We 
also found that the imaginary part of the amplitude of the production of a real photon is 
larger than the one in the case of DIS in a broad range of Q 2 for the reasons as discussed 
above. We also found the same to be true for the full amplitude at small x. We also make 
experimentally testable predictions for the number of real photon events and suggest that 
the number of events are small but not too small such that after improving the statistics 
on existing or soon to be taken data, it would be feasible to access the nondiagonal gluon 
distribution at small x from this clean process. Finally, we demonstrated that measuring 
the asymmetry A at HERA, which is fairly sizable in the kinematics in question, would 
allow one to determine the real part of the DVCS amplitude, in other words gain a first 
experimental insight into nondiagonal parton distributions, despite D < 1. 
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FIG. 4. a) D is plotted versus -/ for x = 1(T 4 , 1(T 3 , 1(T 2 , Q 2 = 12 GeV 2 , B = 5 GeV" 2 and 
y = 0.4. The solid curve is for x = 10" 4 , the dotted one for x = 10 -2 and the dashed one for 
x = 10~ 3 . b) D is plotted versus y for the same x, Q 2 , B and —t = 0.1 GeV 2 
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FIG. 5. a) The asymmetry A is plotted versus 
curve) and x = 1CT 3 (dashed curve) again for Q 2 - 
plotted versus y for the same x,Q 2 ,B and —t = 0. 



t for x = 10 4 (solid curve), x = 10 2 (dotted 
12 GeV 2 , B = 5 GeV~ 2 and y = 0.4. b) A is 
GeV 2 . 
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D plotted versus y for — t = 0. 1 GeV 2 , 
Q 2 = 3.5 GeV 2 . x=10" 4 and 1 0" J 




FIG. 6. a) D is plotted versus -t for x = 1(T 4 and 1(T 3 , Q 2 = 3.5 GeV 2 , B = 8 GeV~ 2 and 
y = 0.4. The solid curve is for x = 1CT 4 , the dashed one for x = 10~ 3 . b) D is plotted versus y for 
the same x,Q 2 ,B and —t = 0.1 GeV 2 
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FIG. 7. a) The asymmetry A is plotted versus —t for x = 10 4 (solid curve), x = 10 3 (dashed 
curve) again for Q 2 = 3.5 GeV 2 , B = 8 GeV -2 and y = 0.4. b) A is plotted versus y for the same 
x,Q 2 ,B and -/ = 0.1 GeV 2 . 
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